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Abstract 

We investigate the simileptonic OZLforbidden annihilation decays B~ — > 
T] {'q,T:^)l~ u for / = /i, e in the perturbative QCD, and carry out a precise 
calculation without any approximation for the one-loop contributions, which 
involves integrals of 4- and 5-point loop functions. Our results show that 
the branching ratios of decays B~ rj l~D, B~ — > rjl~D and B~ t:^1~D, 
turn out to be of orders 10~^, 10~^ and 10~^, respectively, which could be 
observable in the future experiments at the LHC. 
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1 Introduction 



Being the lowest bound state of two heavy quarks(charm and bottom) with open 
(exphcit) flavors, the Be meson provides a unique window into the heavy quark 
physics, which has caused wide experimental and theoretical investigations. After 
the report of the CDF collaboration on the observation of the Be ground state at 
the Fermilab Tevatron people believe that it is possible to accumulate more 
Be meson events in the experiments at the Tevatron Run II|2][3j and the future 
CERN Large Hadron Collider (LHC). At the CERN LHC with the luminosity of 
about C ~ lO^^cm^'^s^^ one can expect around 5 x 10^° Be events per year|3]. 
Unlike symmetric heavy quarkonium (cc and bb bound states), the Be meson is 
composed of heavy quarks with different flavors and it lies below the i?Z)-threshold, 
so its decays via strong and electromagnetic interactions are forbidden. Therefore, 
the investigation of the Be meson decay can offer special information compared to 
symmetric heavy quarkonium. In the framework of the SM its decays can occur 
via three mechanisms: (1) the c-quark decay with the b-quark being a spectator, 
(2) the b-quark decay with the c-quark being a spectator, (3) b-quark and c-quark 
annihilation. The first two mechanisms are expected to contribute about 90% of the 
total width, and the remaining 10% is owed to the annihilation process. 

There is another decay mode which does not belong to the three aforementioned 
types, and it can only occur via the OZI processes. As we know that the Okubo- 
Zweig-Iizuka (OZI) rule|5j plays an important role in the processes which occur via 
strong interaction and in general at the parton level the concerned calculations are 
carried out in the framework of the perturbative QCD (PQCD). Thus careful studies 
on these OZI-forbidden processes can deepen our insight to the perturbative QCD. 
But these processes are very difficult to evaluate, because one not only needs to carry 
out complicated loop calculations at the parton level, but also has to deal with the 
non-perturbative effects of the QCD, which is involved in the hadronization of par- 
tons. Because of lack of solid knowledge on the non-perturbative QCD, for the whole 
calculation one has to adopt concrete models which may contaminate the theoretical 
results. The factorization scheme is just to properly separate calculable perturba- 
tive part from the non-perturbative contributions which must be evaluated either in 
terms of experimental data or using concrete models [6j. In estimating the B decays, 
model-dependent wave-functions are adopted for the light mesons which are the final 
decay products. Among most of the commonly used wave-functions, the light-cone 
wave-function is more favorable, because the finally produced mesons are light and 
so that should be more relativistic. Obviously the errors brought out by using such 
model-dependent wave-function are not accurately estimated, but in some sense they 
are controllable while one uses as much as possible information from available data. 
In order to reduce the theoretical uncertainties of the non-perturbative QCD, the 
meson wave functions should be well theoretically investigated and experimentally 
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tested. To gain more knowledge about the whole picture, one hopes to make the 
part which can be calculated in the framework of the perturbative QCD, as accurate 
as possible and to employ more reasonable model-dependent wave-functions which 
have been tested by fitting data obtained from precise measurements to gain final 
results. Comparing the results with experimental data, one may obtain information 
about both governing mechanisms which are calculated in perturbative framework 
and the wave-functions. In the case, we can also find a trace of new physics. 

The OZI-forbidden process B~ r] l~u was studied by Sugamoto and Yang years 
ago [7j . In their work, an effective Lagrangian was adopted to avoid introducing the 
Be meson wave function, meanwhile they dealt with the light meson by using an 
effective g^g^ rj coupling [7] [Hj, which was obtained in the NRQM(non-relativistic 
quark model) approximation. The valence quark q and anti-quark q in the light 
meson were assumed to possess equal momenta and be on their mass shells, i.e., Pg = 
Pq and = 'TT'q- With such approximation and kinematic assumption, they tactfully 
reduced the complicated Feynman four-point and five-point integral functions P|ITUj 
to be expressed in three-point functions. As for the heavy meson Be, the authors 
neglected the relative momentum and binding energy of heavy quark constituents. 
Then they obtained a reasonable branching ratios of B~ t] Furthermore they 
estimated the branching ratio of B~ 7r°/~z/ based on J/ip decay. 

In this paper, we will investigate the OZI forbiddenjS] B~ rare decay of B^ —>■ 
7] {ri,iT^)l~u in the perturbative QCD without any effective coupling approximation 
and kinematic assumption in one-loop Feynman diagram evaluation. While calcu- 
lating the amplitudes of these decay processes, we have to properly deal with the 
dynamics of bound states. As for the B^ meson, we adopt the way similar to that 
shown in Ref. 0, namely, we ignore the relative momentum of the two heavy con- 
stituents and their binding energy in comparison with their masses. Furthermore, 
it is commonly assumed that the constituents of Be are on mass shell and move 
together with the same velocity for simplicity. In dealing with the final light me- 
son states, we keep an arbitrary relative momentum for the light quarks q and q 
and take the valence quarks (anti-quarks) of the light mesons to be on their mass 
shells. Due to the mass difference between up, down and strange quarks, the decay 
widths of processes B~ —>■ rj (77, 7r*^)/~z/ would receive non-zero contributions. By 
this method, obviously, the advantages of simplifying the loop integral calculations 
vanish. Namely, one cannot reduce the five-point and four-point loop integral func- 
tions into simple expressions involving three-point integral functions. However, with 
the progress in the technique of calculating loop diagrams, we can directly calculate 
the five-point and four-point one-loop integrals numerically in general case without 
any additional approximation and kinematic assumption Then we can obtain 
non-zero rates for the decay modes which obviously violate isospin. 

The paper is organized as follows: In Section 2, we derive the amplitudes. The 
numerical results of the decay rates for the B^ rj {ri,7[^)l^u processes are pre- 
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sented in Section 3, along with all the necessary parameters being listed explicitly. 
Finally, the conclusion is drawn in the last section. 



2 Calculation 



2.1 The amplitude at parton level 

Being OZI-forbidden processes, there is no contribution at the tree-level for the 
decay channel B~ — > Pl~i?, and at one-loop level six diagrams which contribute 
to the decay width, are shown in Fig.l. Here P stands for a pseudoscalar meson 
7] 7r°), and / = /i or e respectively. The six diagrams can be divided into three 
parts as Fig.l (a,d),(b,e) and (c,f). It is clear that the figures in Fig.l(d),(e) and (f) 
can be obtained by exchanging two internal gluons of the corresponding diagrams 
in Fig.l(a),(b) and (c), respectively. 



(a) 



(b) 



(c) 



CFsd' 'CBtJ' ..CMl 

(f) 



(d) 



(e) 



Figure 1: The diagrams at one- loop level for the process — > PI v, where P 
stands for a pseudoscalar meson rj [r], 7r°) 



According to the Feyman rules, the amplitudes at parton level corresponding 
to Fig.l(a), (b) and (c) can be written in the forms as shown in Appendix A. 
By introducing necessary notations(see Appendix B) and carrying out tedious, but 
straightforward derivations, we obtain the explicit expressions of the amplitudes for 
Fig.l(a)-(c). For Fig. 1(a), it is 

iTT^gtg^V.hT^T'' ® T^T" 

KA — ii£ _ 

8[(pi +P2-p-i- Pa? - Ml] [(p2 - P3 - Pa? - m2](27r)4 

X (P3 +Pa- P2)eDl^{u, mq)\Uq{p^)'^^'^5Vq{pA)Vc{p2)lplbUb{pi) 

x[7^(l -75)1/;. (1) 
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The amplitude for Fig. 1(b) is expressed as 



Mb 



8[(pi +P2-P3- PaY - M^] [{p, - p3 - P^y - m2](27r)4 

X(pi - P3 - P4)iDf{u,mq)]Uq{p3h^-f5Vq{p4)Vc{p2hp75Ub{Pl) 



(2) 



and the amphtude corresponding to Fig. 1(c) reads 



8[(pi+P2-P3-P4)2-M2,](27r)4 
xie^„ps[mbm,{X''^''P - z£'^'^"'')E,^(m, m,) - g,^gxe{Y'"'^'^''P + Z^'^^^"") 

X Ef^ {u,mg)]Ug{p3)-f^-f5Vg{p4)Vc{p2) lpl^Ub{pi ) 

xr7^(l -75)z/i. ( 



(3) 



The notations D'^{u, nig), D'^^{u, nig), D'^{u, nig), D^^{u, mq), E"{u, rriq) and E^^^{u, rriq), 
are defined as the one-loop vector/tensor integrals of four- and five-point functions jH] [101, 
integrating over the loop internal momentum k and their explicit expressions can be 
found in Appendix B. 

The definitions of the variables X, Y and Z can be known from the following 
identities, 



The expressions of the contributions from other three diagrams(Fig.l(d)-(f)) are 
similar to that of the first three (for Fig.l(a)-(c)). For convenience, we omit their 
explicit expressions in the text. It is noted that the contribution of Fig. 1(c) involves 
five-point tensor integration functions. We follow the approach in Ref. |9j| to calculate 
directly the five-point scalar and tensor integrals. 




X-'->^^^ - (4) 



'y ^ ^ 'y 'y 




(5) 
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2.2 The hadronic matrix elements 



In above subsection we derive the amplitudes corresponding to the Feynman dia- 
grams shown in Fig.l at the parton level. In order to obtain the decay rates, one 
has to evaluate the hadronic matrix elements. Hadronization happens at the energy 
scale of Aqcd which is the region governed by the non-perturbative QCD. So far, 
there is no any reliable way to evaluate the hadronic matrix elements from any un- 
derlying theory. Instead, to do this job, one needs to invoke concrete models. The 
initial meson Be is composed of only heavy quarks, so that we can suppose the two 
valence quarks b and c in meson to be on mass shell approximately and the com- 
position of Be can be well described by its wavefunction at origin. On the contrary, 
the produced pseudoscalar meson is composed of light quark and antiquark whose 
three-momentum are much larger than Aqcd, so that they are very relativistic, in 
this case, the light-cone wavefunctions seem to be applicable for the hadronization 
of light mesons [m US CSl- 

For the pseudoscalar mesons, the SU(3) flavor wavefunctions are 

r. uu — dd dd + uu + ss dd + uu — 2ss , , 

" = —TT' "° = — 7S — = 7i — ■ 

r] and r]' are mixtures of rjo and rjs, and 

r/ = cos 6*778 — sin 6'?7o, V = smOrjs + cosOrjQ. (7) 

Among the final produced light mesons concerned in the decays of Be — > vr°(r7, rj') + 
W, only rjo is SU(3) singlet. Therefore, the decays Be n^ + W and Be —>■ ris + W are 
related to isospin or SU(3) violation. In general, there are two possibilities which can 
induce isospin or SU(3) violation. The first is photon emission and the second is due 
to mass splitting of u and d quarks. In our case, the violation comes from the quark 
mass splitting in the effective Hamiltonian. There are two sources which are related 
to the quark mass and can contribute to the amplitude. One is from the Feynman 
diagrams where the quark propagators contain quark masses, and another is from 
the higher-twist parts of the wavefunction of meson. But normally we calculate the 
hadron matrix elements with only the leading-twist part of light pseudoscalar wave 
function with valence quarks on-shell, since it can give enough accuracy and simplify 
our calculations. 

The hadronic matrix element part for process B~ Pl~v can be written as 

P{Pp)\qial^q ^ Cr^{u, mg)c-ff3-f5b\Beip^ 
i 

-ifpPpa [ du(j)p{u,n)y2cr'^{u,mg)if^^pg^^, (8) 

JO 
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where P represents 7r°, t] and t]', Ci are the coefficients, the summation is over the 
diagrams shown in last subsection, and in Eq.(jH} we used the expression of the 
leading- twist light-cone distribution amplitudes for pseudoscalar mesons(7r°, rj and 
Tj') with flavor content (gg) written as[T2j 

where the "bar" -notation over u is defined as m = 1 — m, and parameter is the 
renormalization scale of the light-cone operators on the left-hand side. Also the 
light-cone wavefunction of the leading twist is normalized to unity: 



1 

dM0,(M,/i) = l, (P = 7r°,r/,r/'). (10) 



The asymptotic distribution amplitude is defined as the limit in which the renor- 
malization scale goes to infinity. The explicit asymptotic forms of the leading twist 
light-cone wavefunctions of the light pseudoscalar mesons can be different. There 
are several typical leading twist light-cone distribution amplitudes for the light pseu- 
doscalar mesons, so far, one cannot determine which one is the most suitable. In 
our following calculation, we take three different types of the leading twist light-cone 
wavefunctions for light meson, which are frequently adopted |12 | IT U ITH | They 
are expressed as 

(j) p -^{u, jj, oo) = Qu{\ — u), (11) 
0^2(M,/i ^ oo) = 30^2(1 -m)2, (12) 
15 

0,,3(M,/i ^ oo) = y (1 - 2uf[l - (1 - 2uf]. (13) 
Finally, we obtain the hadronic matrix elements of M.ai M-b and A^c as follows 



q=u,d,! 



V Be / 9 q=u,d,s 

xr7^(l -75)i/j, (14) 



7 



^ € q=u,d,s 

x[7^(l -75)z/;, (15) 



xr7^(l -75)z/,. (16) 



Then we can get 

dBr{B- Plu) 



■I-^IV (17) 



dsdt (27r)3 32M3^ 

where is the total amphtude for the process Pl^u at one- loop level and 

the Mandelstam variables s = {pp +piY, t = {pi +PuY- 



3 Numerical results 

In our calculation, no ultraviolet (UV) divergence appears, but there is a mild su- 
perficial infrared(IR) divergence problem. Our method to check the cancellation of 
the IR divergence is standard. Namely, we assign a small gluon mass as a regulator 
and vary it to check if the result is stable. In practical calculations, we set the small 
gluon mass varying in the range between lO""^ MeV to 10"^ MeV and find that 
the result changes only with a negligible small fraction. Therefore, we can trust the 
obtained result which is free of IR problem. In following numerical calculation we 
set the gluon mass being 10~^ MeV . 

The input parameters which we are going to use in the numerical computations 
are taken as follows |7IIIIlII7I[IHlIin] EniEIlES:/^, = SOOMeV", U = ISlMeK, 

= f^, = l57MeV, = 134.9766Me1/, = 5A7.75MeV, = 957.78Me1/, 
rub = 4800MeV^, a,(mBj = 0.20, = ISOOMeK, Mb, = GSOOMeK, Kb = 0.04, 
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tb^ = OAGps, the mixing angle of rj, r/' 6 = —11°, and three possible leading twist 
light-cone wavefunctions of pseudoscalar light meson are given in Eqs. (lllll2ll3|) . 
During the calculation, we keep an arbitrary relative momentum for the light valence 
quark(antiquark) g(g), and moreover, q{q) is on mass shell, while dealing with the 
light meson. 
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rrocess 


rriu 


rrid 




Wr X BR{(pp j 


Wr X BR{(pp,J 


-1 r\f\ T~> T~i { i \ 

Wr X BR{(pp J 




1.5 


4 


/ 

/ 


0.92818 


0.37502 


9.2459 




2 


4 


/ 

/ 


0.50312 


0.21671 


6.2409 




3 


5 


/ 

/ 


0.61314 


0.28907 


5.6131 




3 


7 


/ 

/ 


1.42393 


0.55988 


10.793 




4 


6 


/ 


0.63533 


0.30333 


5.4237 




ITLu 




rris 


10^ X BR{(j)p J 


10^ X BR{(j)p^) 


10^ X BR{(j)p^) 




2 


4 


80 


0.39826 


0.21839 


6.6085 




2 


5 


90 


0.41058 


0.23082 


6.7120 


B- r]l-u 


3 


5 


100 


0.42533 


0.24725 


6.7201 




2 


6 


110 


0.42835 


0.25345 


6.96826 




4 





120 


0.43766 


0.27107 


7.17638 




3 


7 


130 


0.44001 


0.29839 


7.3408 




rriu 


md 


nis 


10^ X BR{(f)pJ 


10^ X BR{(t)pJ 


10^ X BR{(t)p.^) 




2 


4 


80 


0.52462 


0.41698 


2.3586 


B- ri'l'v 


2 


5 


90 


0.52065 


0.40817 


2.3007 


3 


5 


100 


0.51300 


0.40567 


2.2829 




2 


6 


110 


0.50178 


0.39619 


2.2713 




4 


6 


120 


0.49432 


0.39048 


2.2498 




3 


7 


130 


0.49133 


0.38767 


2.2233 



Table 1: The branching ratios of the decays B^ ir^l^u, B^ — > rjl^u and B^ — >■ 
1] l~i> in the rest frame of Be are hsted, and the three columns correspond to the 
three different leading twist light-cone wavefunctions of the produced pseudoscalar 
light mesons {n^, t], rj'). The masses for the light quarks(u, d, s) are in MeV. 

We carry out the integration of the scalar and tensor four- and five-point integral 
functions precisely. We adopted the FF package [SB] in the calculation of two-, three- 
and four-point integral functions, while the implementations of the scalar and the 
tensor five-point integrals are done exactly by using the Fortran programs as we 
used in our previous works on e~^e~ — >■ ttH^ and e^e~ — > Z^H^H^ processes PH 123 
by using the approach presented in Ref.[9 . 

In order to model the light hadronic effects, we use three different types of wave- 
functions for light mesons in our calculation. We present the theoretical predictions 
on the decay widths of B~ — > tt^I^u, B^ — > rjl'i/ and B^ ^ r] l~i> in the rest frame 
of Be, corresponding to the three different leading twist light-cone wavefunctions of 
light mesons respectively in Table 1. 
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4 Conclusion 



In this work, we studied the OZI forbidden quark-level sub-processes of — >■ 
1] {ri,TT^)l~'iJ in the framework of the perturbative QCD. For B~ 7r°/~z/ process, 
within reasonable ranges of the masses of u and d quarks, the branching ratio with 
the distribution amplitude (pp^ can be about one order larger than other two dis- 
tribution amplitudes {(pp^.^ip^) and can reach the order of 10~^ with specific light 
quark masses. For B~ — > rjl'v process, the branching ratios turn out to be of the 
order 10~^ for distribution amplitudes and 0^2, and about 10~^ for distribu- 
tion amplitude (pp^. For B^ ri'l^u decay, the (j)p^ can give the branch ratio 
values of the order 10~^. These results seem to be at the reach of future experi- 
ments at the LHC where it is expected to produce 5 x 10^° Be meson events per 
year with ^/s = 14 TeV and luminosity C = 10^^cm~^s~^|5. We conclude that the 
decays B^ r] {ri,Tc'^)l~u could be investigated experimentally at the LHC, and 
the study on these decays can deepen our understanding on both perturbative and 
non-perturbative QCD. 
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Appendix A 

The amplitudes at parton level corresponding to Fig.l (a), (b) and (c) are ex- 
plicitly expressed as 



2^/2 + ky - ky {pi + P2 - P3 - pa)^ - 



aU rparpb ^ rnbrna ^ I a hv 



W 

8 [pi +P2-P3- PaY - J (27r)4 

xvq{pi)vc{p2)l'^{,i>A + "^c)7"(l^3 + 'mc)l^(l - 75)Ufe(pi)r7^(l 

1 



- m2)(p3 + kY{pi - ky{{p2 -P4 + ky - ml){{p2 - Ps - PaY - m^) ' 
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X 



^cb 

-f-rug - 
-igsTV)- i-igsT^r)uM 



i)^- i>i- rrib fx - i>z - - 



1 

- m2)(p3 + /i;)2(p4 - - ps - k)^ - m^)((pi - P3 - PaT - ml) ' 



p^%fe)(-»g.r''7») _^ ^ (-ig.r''7,3)f,(P4)iic(P2)(-»g.rV) 

xZ— 17^^(1 - 75)1/; X 



2^/2 {P3 + kf (p4 - ^? {Pi +P2-P3- Pa? - M^' 



X'yg(]94)wc(p2)7 (A -|^2 + ?"c)7^(l - 75) (A -'h-^ + mb)-f''ub{pi)hfj,{l - 

X \ (20) 

- m2)(p3 + /c)2(p4 - - p^ - k? - mg)((p4 - P2 - k? - ml) ' 



Appendix B 

The notations of the four-point and five-point one-loop functions in our text are 
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defined as 



If k'^ 
^« = ^ y "^'^[F - m2](p3 + k){p, - ky[{p2 - P4 + - m2] ' ^^^^ 
m) = — f d'^k k-{p,-p2-kf 

1 /" k'^ 

DI{U, m,) = — y ^'^[^2 _ ^2](p3 + _ fc)2[(p3 - - ^)2 _ ^2] ' (23) 
j^ai( ^__}_ [ M, k'^jpi -P3- k)^ 

' [k^-m^Mp, + k){p,-kn{p,-p,-kf-mlY ^'^^ 



k" 



[A;2 - m2](p3 + k){pi - ky[{pi -p2 - k^ - mf\[{pi - ps - k) - ml] ' 

(25) 

k''{p4-P2-k)^ipi-p3-ky 



"^q](P3 + k){p4 - ky[{p4 -P2- kf - m2][(pi - P3 - A;) - mg] ' 

(26) 



with 



777,6 ?T7c . 

Pi = IT^— PSo P2 = ITj—P^-' ^ ^4 = Pp(l - «)■ 
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